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Abstract 

Taking into account the A baryon distribution amplitudes and the most general 
form of the interpolating current of the A^, the semileptonic A;, — )• Ai^i^ transition 
is investigated in the framework of the light cone QCD sum rules. Sum rules for 
all twelve form factors responsible for the A^ — )• Ai^i~ decay are constructed. The 
obtained results for the form factors are used to compute the branching fraction. A 
comparison of the obtained results with the existing predictions of the heavy quark 
effective theory is presented. The results of the branching ratio shows the detectability 
of this channel at the LHCb in the near future is quite high. 
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1 Introduction 



Experimentally, the detection and isolation of the heavy baryons is simple comparing to the 
light systems since having the heavy quark makes their beam narrow. In the recent years, 
considerable experimental progress has been made in the identification and spectroscopy 
of the heavy baryons containing a heavy bottom or charm quark [1-8]. These evidences 
can be considered as a good signal to search also the decay channels of the heavy baryons 
hke, Afe — > kt^i" at LHCb. This rare channel, induced by the flavor changing neutral 
currents (FCNC) of fe — > s transition, serves testing ground for the standard model at loop 
level and is very sensitive to the new physics effects [9], such as supcrsymmctric particles 
[10], light dark matter [11] and also fourth generation of the quarks and extra dimensions, 
etc. Moreover, this channel can be inspected as a useful tool in exact determination of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, and Vts-, CP and T violations, 
polarization asymmetries. 

Theoretically, there are some works devoted to the analysis of the heavy baryon decays, 
where practically in all of them the predictions of the heavy quark effective theory (HQET) 
for form factors have been used. Transition form factors of the A;, — )■ Ac and Ac ^ A decays 
have been studied in three points QCD sum rules in [12], the A;, — )■ plv transition form 
factors have also been calculated via three point QCD sum rules in the context of the heavy 
quark effective theory (HQET) in [13] and in the framework of the SU(3) symmetry and 
HQET in [14]. In the present work, using the most general form of the interpolating current 
for the A;, and also the distribution amplitudes of A baryon, all form factors related to the 
electroweak penguin and weak box diagrams describing the A;, — > Ki'^t' are calculated in 
the frame work of the light cone QCD sum rules in full theory. The obtained results for 
the form factors are used to estimate the decay rate and branching ratio. Regard that this 
transition has been investigated in [15] and [16] also in the context of the HQET but the 
same frame work using the distribution amplitudes of the A and A;,, respectively. Moreover, 
form factors, branching ratio and di-lepton forward-backward asymmetries are studied in 
[17-19] also within the context of the HQET. In [20-22], E^^c and A^^c to nucleon transitions 
are also evaluated using the nucleon wave functions in light cone QCD sum rules approach. 

The plan of the paper is as follows: in section II, the light cone QCD sum rules for the 
form factors are obtained using the A DA's. The HQET relations among all form factors 
are also discussed in this section. Section III is dedicated to the numerical analysis of the 
sum rules for the form factors as well as numerical results of the decay rate and branching 
ratio. 

2 Theoretical Framework 

The Afe — )■ A£+£^ channel proceeds via FCNC 6 — )■ s transition at quark level. The effective 
Hamiltonian describing the electroweak penguin and weak box diagrams related to this 
transition can be written as : 
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2mb C-!— sia^^q^il + 75)6/7''^ 



(1) 



To find the amplitude, we need to sandwich this effective Hamiltonian between the initial 
and final baryon states, i.e., (Afe(p + g)|'He//|A(p)). From Eq. (1) we see that in calculation 
of the Af, — > Ai~^i~ decay amplitude, the matrix elements, (Aft(p + g)| 67^^(1 — 75)s|A(p)) and 
{Ai,{p + q)\bia^u(l'^{l + 'y5)s\A{p)) arc appeared. These matrix elements can be parametrized 
in terms of the twelve form factors, fi, Qi, ff and gj in the following manner: 



(A(p) I ^7^.(1 - l,)h I A,(p + q)) = ua{p) [^MQ') + ^^Mv9V2(g') + q'MQ') 

- lnl59iiQ^) - i(ytJ,uim''g2{Q'^) - q^l^g^iQ"^)] UA,ip + q) , (2) 



and 



+ wqUQ^) + i(y^ivi^q"gl{Q^) + q'"i^gl{Q'')\ ua,{p + g) , (3) 

where = — g^. For calculation of these form factors we use the QCD sum rules approach. 
To obtain the sum rules for the form factors in this approach, the following correlation 
functions, the main objects in this approach, are considered: 



nj(p,g) 



/ 



d'xe-'^^{0 I T{J^H0),&(^)7m(1 -75)s(^))} I Hp)) , 



(4) 



where, p represents the A's momentum and q is the transferred momentum and the J^'' 
is interpolating current of A^. The most general form of the interpolating current of A^ 
baryon can be written as: 



J^'ix) 



V6 



{qr {x)Cq',{x))^,b'={x) + /3{qr {x)C^,q',{xW{x) 



+ {qr {x)Cb\x))^,q^,{x) + ^{qr {x)C^,b\xM{x) 



(5) 



where qi and q2 are the u and d quarks, respectively, a, b and c are color indexes and C is the 
charge conjugation operator. The ^ is an arbitrary parameter with ^ — —1 corresponding 
to the loffc current. 

In order to obtain the sum rules for the transition form factors, we will calculate the 
aforementioned correlation functions in two different ways, namely, physical (phenomeno- 
logical) and theoretical (QCD) sides and equate these two representations isolating the 
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ground state through the dispersion relation. Finally, to suppress the contribution of the 
higher states and continuum, we will apply the Borel transformation and continuum sub- 
traction to both sides of the correlation function and impose the quark hadron duality 
assumption. 

The first step is to calculate the physical side of the correlation functions. Saturating the 
correlation functions with complete set of the intermediate states with the same quantum 
numbers as the initial state, for the physical part of the correlation function we obtain, 

n^(p _ y (0 I -^^"(o) I Hp + s)){Hp + I ^7m(i - ts)^ I a(p)) ^ .g. 



(0 I J-^''(0) I A.Q; + q. s)){A,{p + q. ,s) | bia,,,jr{l + - ,),s | A(p)) 

^A. - (P + 



+ ••• , (7) 



where the dots • • • represent the contribution of the higher states and continuum. The 
vacuum to the baryon matrix element of the interpolating current, (0 | J^''(0) | Ah{p + q, s)) 
is written in terms of the residue, Aa^ as: 



(0 I J^'(O) I Ab{p + q, s)) = XaqUaJp + q,s 



(8) 



Putting Eqs. (2), (3) and (8) in Eqs. (6) and (7) and performing summation over spins of 
the Aft baryon using 



^UA^{P + Q, s)uA^{p + q, s) ^+mA, , 



(9) 



we get the following expressions for the correlation functions 

^+ ^ + rriA, 



n>,g) = A 



^'ml - {p + qf 



n,(p,g)-AA. ^, -(p + g)^ 



+ iK.i5gI + icf^iuibq^'gl - q^ibgl'^UA{p) . 



(10) 



(11) 



Using the equation of motion and Eqs. (10) and (11), we get the following final expressions 
for the phenomenological sides of the correlation functions: 



K^p.q) 



A 



At 



:{Vl{Q^)p^. + V2{Q^)p, i+ [/2(Q') + /3(Q')]g^ i 



- {p + qY 

- 2yi(g2)p^75 + 2y2(g')PM ^75 + U(g') + g^{Q')\q^. ^75 



+ other structures |ma(p) , 



(12) 



A/ 



- (p + q)^ 



+ other structures I tiA(p) . 



(13) 



To compute the form factors or their combinations, /i, /2, /2 + /a, 5'2 and g2 + (73, 
we will choose the independent structures p^, p^^, g^^, ]?^75, Pi^^jb, and 5^1x^75 from Eq. 
(12), respectively. The same structures are selected to calculate the form factors or their 
combinations labeled by T in the second correlation function, Eq. (13). 

The next step is to calculate the correlation functions from QCD side in deep Euclidean 
region where {p + g)^ ^ 0. For this aim, we expend the time ordering products of the 
interpolating current of the A(, and transition currents in the correlation functions (see 
Eq. (4)) near the light cone, ~ via operator product expansion, where the short 
and long distance effects are separated. The former is calculated using QCD perturbation 
theory, whereas the latter are parameterized in terms of the A DA's. Mathematically, this 
is equivalent to contract out all quark pairs in the time ordering product of the J^'' and 
transition currents via the Wick's theorem. As a result of this procedure, we obtain the 
following representations of the correlation functions in QCD side: 



2(C)^^(75)p^ + (C)^^(75)p0 + (C);30(75) 



2(C75),</,(/)p/3 + {C'y,U{I)p^ + (C75)/3<a(/),p] } [7^(1 - 75) 



X 56(-x)/3.(0|<(0)s^(x)d^(0)|A(p)) , 



(14) 



+ 

X Shi 



{ [2(C)^0(75)p/3 + iC)r^/3{l5)p4> + iC)p^{-f5)r,p 



2{C-i^\^{I)pp + {C-i^U{I)p^ + (C75);3^(/),p] } 

-x)pMK{^)s'e{x)dimm) , 

The heavy quark propagator, Sb{x) is calculated in [23]: 



iq^a^^il + 75) 



(15) 



Sb{x) = Sl'^'^x) - ig. 



d^k 







^G^''{vx)a^, 



{rriQ - k'^y 



rrin — k^ 



(16) 



where. 



gfree 



4^ 7^ ' 



i.Tl'^x'^ 



K2{mbV-x'^) , 



(17) 
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and Ki are the Bessel functions. The terms proportional to the gluon field strength are 
contributed mainly to the four and five particle distribution functions [23-27] and expected 
to be very small in our case, hence when doing calculations, these terms are ignored. The 
matrix element e"^'^(0|it^(0)(ig(0)s^(x)|A(p)) appearing in Eqs. (14,15) represents the A's 
wave functions, which are calculated in [27] and we list them out for the completeness of 
this paper in the Appendix. Using the A wave functions and the expression of the heavy 
quark propagator, and after performing the Fourier transformation, the final expressions 
of the correlation functions for both vertexes are found in terms of the A DA's in QCD or 
theoretical side. 

In order to obtain the sum rules for the form factors, /i, f^, fs, gi, g2, Qs, ff, /J, /J, 
gj, gj and g^, we equate the coefficients of the corresponding structures from both sides of 
the correlation functions through the dispersion relations and apply Borcl transformation 
with respect to {p + q)^ to suppress the contribution of the higher states and continuum. 
The expressions for the sum rules of the form factors are very lengthy, so we will give only 
extrapolation formulas to explore their dependency on the transferred momentum squared 

The explicit expressions of the sum rules for the form factors depict that to calculate 
the values of the form factors, we need also the expression of the residue, Aa(,- This residue 
is calculated in [28]. 

Few words about the relations among the form factors in HQET are in order. In HQET, 
the number of independent form factors is reduced to two, Fi and F2, so the transition 
matrix element can be parameterized in terms of these two form factors as [28,29]: 

(A(p) I sFb I At,{p + q)) = UA{pm{Q')+ j^F2{Q')]ruj,,{p + q), (18) 

Here, F refers to any Dirac matrices and ^ = (^+ i[)/m\^. Comparing this matrix element 
with definitions of the form factors in Eqs. (2) and (3), the following relations among the 
form factors are obtained (see also [30,31]): 

h-gi-fI-gI-F, + ^F,, 

F2 

f2 = g2 = fs = gs = — , 

rT T ■P'2 2 

fl = 9i ^ Q , 

fl = -{rriA^ - ^a) , 

F2 

gj = ("^A(, + mA) . (19) 



3 Numerical Analysis 

This section is devoted to the numerical analysis of the form factors, their extrapolation 
in terms of the momentum transferred square and calculation of the total decay rate and 
branching ratio for rare A^ — > A£'^£~ transition in QCD. 



5 



Some input parameters used in the numerical calculations are: {uu){l GeV) = {dd){l GeV) — 
-(0.243±0.01)3 GeV'^, mg(l GeV) = (0.8±0.2) GeV^ [32], ttia = (1115.683±0.006) MeV, 
niA^ = (5620.2 ± 1.6) MeV and = (4.7 ± 0.1) GeV. Sum rules for the form factors de- 
pict that the A DA's are the main input parameters (see the Appendix). They contain 4 
independent parameters which are given as [27]: 

/a = (6.0 ± 0.3) X 10-3 GeV^ , Ai = (1.0 ± 0.3) x 10"^ GeV^ , 

IA2I = (0.83 ± 0.05) X 10-^ GeV^ , [Agl = (0.83 ± 0.05) x 10"^ GeV^ . (20) 

It is well known that, the Wilson coefficient Gg^^ receives long distance contributions 
from J/ip family, in addition to short distance contributions. In the present work, we do 
not take into account the long distance effects. From the explicit expressions for the form 
factors, it is clear that they depend on three auxiliary parameters, continuum threshold 
So, Borel mass parameter and the parameter (3 entering the most general form of the 
interpolating current of the A^. The form factors should be independent of these auxiliary 
parameters. Therefore, we look for working regions for these parameters, where the form 
factors are practically independent of them. To determine the working region for the Borel 
mass parameter the procedure is as follows: the lower limit is obtained requiring that the 
higher states and continuum contributions constitute a small percent of the total dispersion 
integral. The upper limit of is chosen demanding that the series of the light cone 
expansion with increasing twist should be convergent. As a result, the common working 
region of is found to be 15 GeV'^ < < 30 GeV'^. As an example, we present the 
dependence of the form factor /i on Borel mass parameter, at two fixed values of 
in Fig. 1. From this figure it follows that the form factor fi exhibits a good stability with 
respect to the variations of M^. The continuum threshold Sq is correlated to the first exited 
state with quantum numbers of the interpolating current of the A;, and is not completely 
arbitrary. Numerical analysis leads to the interval, (mA^, + 0.3)^ < sq < (mA^ + 0.5)^, where 
the form factors weakly depend on the continuum threshold. In order to attain the working 
region for the parameter, /3, we look for the variation of the form factors with respect to 
cos^, where f3 = tanO. After performing numerical calculations, we obtained that in the 
interval —0.6 < cos^ < 0.3 all form factors weakly depend on j3. As an example, we 
show the dependence of the form factor, /i on cos 9 at two fixed values of the and at 
= 22 GeV^ in Fig. 2 . From this figure indeed we see that in the aforementioned region 
of cos^^, the form factor fi weakly depends on /3. 

The analysis of the sum rules, as has already been explained above, is based on, so 
called, the standard procedure, i.e., the continuum threshold Sq is independent of and 
q^. However, in [33], instead of the standard procedure, namely, independence of the sq 
from Af^ and g^, it is assumed that the continuum threshold depends on Af^ and q^ and 
this leads to large realistic errors. Following [33], in the present work the systematic error 
is taken to be around 15%. 

In calculating the branching ratio of the A^ — > Kt^l' decay, the dependence of the form 
factors fi{q^), gi{q^), ffiq^), and gj{q^) on q^ in the physical region Ami < <f < {mA^-rriAY 
are needed. But unfortunately, sum rules predictions for the form factors are not reliable 
in the whole physical region. Therefore, in order to obtain the dependence of the form 
factors from sum rules we consider a range of q^ where the correlation function can reliably 
be calculated. For this aim we choose a region which is approximately 1 GeV below the 
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QCD sum rules 




a 


b 




fl 


-0.046 


0.368 


39.10 




0.0046 


-0.017 


26.37 


h 


0.006 


-0.021 


22.99 


91 


-0.220 


0.538 


48.70 


92 


0.005 


-0.018 


26.93 


93 


0.035 


-0.050 


24.26 


fl 


-0.131 


0.426 


45.70 


fl 


-0.046 


0.102 


28.31 


9l 


-0.369 


0.664 


59.37 


9l 


-0.026 


-0.075 


23.73 



Table 1: Parameters appearing in the fit function of the form factors, /i, /2, /s, g^^ gs, 
/J) fl: 92 9I f^ll theory for — > M^l~ . In this Table only central values of the 
parameters are presented. 



perturbativc cut, i.e., up to ~ 12 GeV"^. To be able to extend the results for the form 
factors to the whole physical region, we look for a parameterization of the form factors in 
such a way that, in the region 4m^ < < 12 GeV'^ this parameterization coincides with 
the sum rules predictions. 

The next step is to present the dependency of the form factors. Our numerical 
calculations show that the best parameterization for the dependence of the form factors /i, 
/2, /s, 91, 92, gs: fl: fl, gj and gj on q^ is as follows: 

fiiQ'mq')] = , " , , + ^-r—, , (21) 



1 - 



1 - 



m 



fit 



where the fit parameters a, b and in full theory are given in Table 1. On the other 



hand, we find that the best fit for the form factors f^ and g-f is of the following form. 



1 



m 



fit 



1 - 



m 



2 ■ 



(22) 



The results for the parameters c, rnj^^ and m^f^ are presented in Table 2. In extraction of the 
values of the fit parameters presented in both Tables 1 and 2, the values of the contimmm 
threshold, sq = 35 GeV^, Borel mass parameter, = 22 GeV^ and cos 6* = 0.2 have been 
used. 
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QCD sum rules 




c 


'2 


"2 


fl 

9l 


-1.191 
-0.653 


23.81 
24.15 


59.96 
48.52 



Table 2: Parameters appearing in the fit function of the form factors and gj in full 
theory for Ai+i~. 

The values of form factors at = are also presented in Table 3. In this table we also 
present the numerical results obtained from HQET, using the values for the form factors 
Fi(0) = 0.462 and F2 = -0.077 predicted in [17], and relations in Eq. (19) at HQET 
limit. The errors in the values of the form factors at = are due to the uncertainties 
coming from M^, sq, the parameter (3, errors in the input parameters, as well as from the 
systematic errors. From this Table we see that, the predictions of the HQET on the form 
factors are changed more than 40% for the form factors /i(0), gi{0), /J(0), and gj{0), while 
the results of both approaches are very close to each other for the remaining form factors. 





Present work 


HQET ([13]) 


/i(0) 


0.322 ±0.112 


0.446 


/2(0) 


-0.011 ±0.004 


-0.013 


/3(0) 


-0.015 ± 0.005 


-0.013 


^i(O) 


0.318 ±0.110 


0.446 


^2(0) 


-0.013 ± 0.004 


-0.013 


^?3(0) 


-0.014 ±0.005 


-0.013 


/r(o) 


0±0.0 


0.0 




0.295 ±0.105 


0.446 


flio) 


0.056 ± 0.018 


0.061 




0±0.0 


0.0 




0.294 ±0.105 


0.446 


gm 


-0.101 ±0.035 


-0.092 



Table 3: The values of the form factors at = for M'^l . 

The final task is to calculate the total decay rate of the A;, — )■ kt^i" transition in the 



8 



whole physical region, 4mj < < {rn\^ — mjCf'. The differential decay rate is obtained as: 



dV 
ds 



-8192^1 WA 



e(s) + -A(s) 



(23) 



where s = /m\ , r = m\/m\ , A = A(l, r, s) = — 2r — 2s — 2rs, Gp = 1.17x 10" 



4m? 



GeV ^ is the Fermi coupling constant and v = y 1 — is the lepton velocity. For the 

element of the CKM matrix | Vj^V^* |= 0.041 has been used [34]. The functions 0(s) and 
A(s) are given as: 

@{s) = 32m^,m\^s{l + r - s) (lAsl' + l^sl') 
+ 6Amlml^{l - r - s) Re[D*^3 + ^3^*] 
+ 64m^^v^(6m2 - m^^s)Re[D*Ei] 

+ 64mjml^^/^(2mA,sRe[D;Es] + (1 - r + s)Re[L'*L>3 + ElEs]j 

+ 32mi^(2m| + mi^s)|(l - r + s)mA,^^Re[>l^>l2 + 5*^2] 

- mA,(l - r - s) Re[^t52 + A^Ei] - 2^? (^Rc[A*5i] + m^^s Re[A;52]) } 

+ 8mi^{4mKl + r - s) + m^^ [(1 - rf - s'] } + \B,f) 



+ 8m^^|4m^ A + (l + r-s)s + m^^s (1 



\2 2 

r — s 



8m 



,i^{4m,2(l + r - s) - mi^ [(1 - rf - s'] } {\D^f + \E,f) 
+ 8m%mp\^ - SmA^sy/¥Re[D*E2] + 4(1 - r + s)^^ Re[D*r>2 + E*E2] 
- 4(1 - r - s) Re[DlE2 + + m^^ [(1 - r)^ - s^j (|i^2|2 + |^2|2) | ^ (34) 



where 



A (s) = -8miyA (lAl' + \B^f + \D^f + \E,f) 
+ 8mi^st.2A(|A2|' + |S2|' + |L>2|' + |^2|') , 



^1 = ^ (A^ + al) {-2m,Cr) + (/i - ^1) a 

A2 = (1 ^ 2) , 
yl3 = (1 ^ 3) , 

5i = Ai {gi -gi; g{ -g^) , 
B2 = fii (1 ^ 2) , 
B, = fii (1 ^ 3) , 
Di = (/i - gi) Cio , 
D2 = (1 ^ 2) , 



eff 



(25) 



(26) 
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Ds = (1 ^ 3) , 

El = Di {gi -gi) , 
E2 = ^1 (1 ^ 2) , 

^3 = ^1 (1 ^ 3) . (27) 

Integrating the differential decay rate on s in the whole physical region A.ml/m\^ < s < 
(1 — \/r)^ and using the hfe time of the A5 baryon, t\^^ = 1.383 x 10~^^ s [34], we obtain 
the results for the branching ratio which are presented in Table 4. 
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HQET([19]) 


Br{Kb Ae+e-) 


(4.6 ±1.6) X 10-^ 


(2.23 ^ 3.34) X 10-6 


5r(Ab Afi+fjr) 


(4.0 ±1.2) X 10-6 


(2.08^3.19) X 10-6 


Br{Ab Kt+t-) 


(0.8 ± 0.3) X 10-6 


(0.179 ^ 0.276) X 10-6 



Table 4: Values of the Branching ratio for A^ — )■ A in full theory and HQET for 

different leptons . 

In this Table we also present the values of the branching ratio obtained in HQET [19]. 
Comparing the results of both approaches, we see that our predictions on the branching 
ratios for the A;, Ae+e-, A;, A/i'^ij,' channels are larger approximately by a factor 
of two than the ones predicted by the HQET, while for the A^ — > At'^t~ channel our 
prediction is four times larger than the result of the HQET. Since 10^° 10^^ pairs are 
expected to be produced per year at LHCb [35], the results presented in Table-4 show that 
detectability of A;, — )■ Ai'^i^ {£ = e, fi, r) decays in this machine is quite high. 

In conclusion, we calculate all twelve form factors responsible for the A^ — )■ A£+£^ decay 
within light cone sum rules. It is obtained the maximum difference between our results and 
HQET predictions on the form factors is about 40%. Using the parametrization for the 
form factors, the branching ratio of the A^ — )■ Ai~^i~ decay is estimated, and the result we 
obtain allows us to conclude that the delectability of this decay at LHCb is quite high. 
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Appendix 



In this Appendix, the general decomposition of the matrix element, e"'"^(0|M^(0)sg(a;)(i^(0) | A(p)) 
entering Eqs. (14,15) as well as the A DA's are given [27]: 

4(0|e"'X(«ia^)4(«2x)rf;(a3a:)|A(p)) 



2 2 

+ nmA(75C)a^A^ + V2ml{^,CU{:f:K)^ + (Vi + '^V^){i)CU{^,K)^ 
+ V2mK{j)C)^p{ic-ir,K)^ + V3mA(7/,C)„^(7'^75A)^ + V4mi(;^C)a^(75A)^ 

2 2 

+ ^p^^f )(|^75C)„/3A^ + A"^a(j^75C)„/3(M)7 + '43mA(7^75C),^(7^A)^ 
+ AmA(ji^75C)a/3A7 + A"^i(7/^75C)a/3(«o-''''x^A)^ + A&m\{^^zC)^p{^h)^ 



+ (71 + 



2^2 



^7^^)(p'^ia^,C)„;3(7^75A)^ + r^rriKixi'p^ia^.CUii.^^y 



+ 7imA(o-^,.C)a/3(o-^''75A)^ + 7imA(p''(T^,.C)a/3(a^^Xp75A)^ 



(Al) 



The calligraphic functions in the above expression have not definite twists but they can be 
written in terms of the Lambda distribution amplitudes (DA's) with definite and increasing 
twists via the scalar product px and the parameters a^, i = 1, 2, 3. The explicit expressions 
for scalar, pseudo-scalar, vector, axial vector and tensor DA's for Lambda are given in 
Tables 5, 6, 7, 8 and 9, respectively. 



2pxS2 = Si - S2 



Table 5: Relations between the calligraphic functions and Lambda scalar DA's. 



Vi = Pi 



2pxV2 



Table 6: Relations between the calligraphic functions and Lambda pseudo-scalar DA's. 
Every distribution amplitude F{aipx)= Si, Pj, V^, Ai, Tj can be represented as: 



F{aipx) = J dxidx2dx^6{xi +X2 + X2,- i)e**'^^'^'''*Fl 



Xi 



(A2) 



where, Xi with i = 1, 2 and 3 are longitudinal momentum fractions carried by the partici- 
pating quarks. 
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Vi = Vi 
2pxV2 = Vi-V2-V3 



2V, = V. 



4pxV4 = -2Vi + V3 + V4 + 2V5 



4pxV5 = V4-V5 



AipxfVe = -Vi + 1/2 + 14 + 14 + 14-14 



Table 7: Relations between the calligraphic functions and Lambda vector DA's. 



Ai = Ai 



2pxA2 



-Ai + - A3 



2As = A. 



ApxAi = -2Ai - A3- Ai + 2A5 



ApxA^ = A3 - A4 



4.{pxyA& = Ai - A2 + A3 + A4 - A5 + A 



Table 8: Relations between the caUigraphic functions and Lambda axial vector DA's. 



The explicit expressions for the A DA's up to twist 6 are given as: twist- 3 DA's: 

Ai{Xi) = -120X1X2X303 , 



V^{x,) = , 
Ti{xi) = . 

twist-4 DA's: 

5i(x,) = 6X3(1 -X3)(e4+el°), 

V2{Xi) = , 

14 (Xi) = 12 (Xi - X2)X3V'4 , 

T2{Xi) = , 

T,{Xi) = -6(xi-X2)x3(e4+el°) 

Ttwist-5 DA's: 

3 



<S'2(Xi) 

V^{x,) 
T^ixi) 



^{xi + X2){e, + u') , 

3(x2 - xOV's , 
0, 



Tsixi) = 
and twist- 6 DA's 



^{x, - X2M - ^',') 



V(i{x,) = , 

nixA = . 



Pi{X:i) 

A2{Xi) 

Asixi) 
Tsixi) 



P2{Xi) - 

Aiixi) 
A^ixi) 

n{xi) 



Aeixi) 



6(i-.x3)(^4°-e;°) 

: -24X1X201 , 
: -12X3(1 -Xs)^. 



4 ) 



6(x2-Xi)x3(-e4° + e?) 



-3(1 - X3)V^° , 

-6a;305 > 
0, 



(A.3) 



(A.4) 



(A.5) 



(A6) 
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Tl=Ty 




2pxT2 = Ti + T2 


-2T3 






2px% = Ti - T2 


- 2T7 




2px% = -Ti + T 


+ 2n 


4:{pxfr6 = 


-- 2T2 - 2T3 - 2T4 - 


f 2T5 + 2T7 + 2T8 




ApxTj = T7 — 


Ts 


4:{pxfrs 


= -Ti + T2 + T5 - 


- Tg + 2T7 + 2Ts 



Table 9: Relations between the calligraphic functions and Lambda tensor DA's. 

The following functions are encountered to the above amplitudes and they can be defined 
in terms of the 4 independent parameters, namely f\, Ai, A2 and A3: 

4 = <l>6=-fA, 0^ = 0^ = -^(/a + Ai) , 

^4 = ^5 = lih- ^1) ' ^4 = ^5° = A2 + A3 , 

= = A3 - A2 . (^.7) 
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Figure 1: The dependence of form factor, /i on Borel mass parameter at two fixed values 
of the gf^, and at Sq = 35 GeV'^ and /3 — 5. 
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Figure 2: The dependence of form factor, /i on cos^ parameter at two fixed values of the 
q^, and at Sq = 35 GeV^ and M| = 22 GeV^. 
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